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ABSTRACT 



Fermilab experiment E791, measuring charmed hadron production in 7r~A inter- 
actions at 500 GeV with high statistics, has observed a strong asymmetry between the 
hadroproduction cross sections for leading D mesons which contain projectile valence 
quarks and the nonleading charmed mesons without projectile valence quarks. Such 
correlations of the charge of the D meson with the quantum numbers of the beam 
hadron explicitly contradict the factorization theorem in perturbative QCD which 
predicts that heavy quarks hadronize through a jet fragmentation function that is in- 
dependent of the initial state. The E791 experiment also measures Ac/ Ac and Dg/Dg 
production asymmetries as well as asymmetries in DD pair production. We exam- 
ine these asymmetries and the fractional longitudinal momentum, xp, distributions 
for single and pairs of charmed hadrous within a two- component model combining 
leading-twist gg and qq fusion subprocesses with charm production from intrinsic 
heavy quark Fock states. A key feature of this analysis is intrinsic charm coalescence, 
the process by which a charmed quark in the projectile's Fock state wavefunction 
forms charmed hadrons by combining with valence quarks of similar rapidities. 
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1. Introduction 



The E791 experiment, studying 500 GeV tt'A interactions with carbon and plat- 
inum targets, employs an open geometry spectrometer with a very open trigger and 
a fast data acquisition system to record the world's largest sample of hadroproduced 
charm [^. This large data set allows detailed investigations of charmed hadron pro- 
duction including Ac{udc), Ac, Ds{cs), Ds, and DD pairs. 

One of the most striking features of charm hadroproduction is the leading particle 
effect: the strong correlation between the quantum number of the incident hadron 
and the quantum numbers of the final-state charmed hadron. For example, more D~ 
than are produced at large xp in ii'A D^X |^, ^ ^, |^. There is also 
evidence of leading particle correlations in Ac |P, |^, and Kbiudh) production in 
pp collisions and Sc('usc) production in hyperon-nucleon interactions [1^, |lT|. Such 
correlations are remarkable because they explicitly contradict the factorization theo- 
rem in perturbative QCD which predicts that heavy quarks hadronize through a jet 
fragmentation function that is independent of the initial state. 

Leading charm production can be quantified by studies of the production asym- 
metries between leading and nonleading charm production. In Ti~p interactions, both 
D~{cd) and D^{cu), which share valence quarks with the n^{ud), are "leading" while 
D~^{cd) and D^{cu), which do not, are "nonleading" at xp > 0. The harder leading 
D distributions suggest that hadronization at large xp involves the recombination of 
the charmed or anticharmed quarks with the projectile spectator valence quarks. The 
D^/D^ asymmetry is defined as 

_ da{D~) - da{D+) 
^'''"'^ - da{D-) + da{D+) ' 

The measured asymmetry increases from nearly zero at small xp to Ad-/d+ ~ 0.5 
around xp = 0.65 0, ^, indicating that the leading charm asymmetry is primarily 
localized at large xp. Thus the asymmetry Ad-/d+ refiects the physics of only a 
small fraction of the total cross section. The neutral D's were not used in the 
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analysis since these states can also be produced indirectly by the decay of nonleading 
mesons. 



In a recent paper [|12|, we discussed a QCD mechanism which produces a strong 
asymmetry between leading and nonleading charm at large xp- A key feature of 
this analysis is coalescence, the process by which a produced charmed quark forms 
charmed hadrons by combining with quarks of similar rapidities. In leading-twist 
QCD, heavy quarks are produced by the fusion subprocesses gg QQ and qq QQ. 



The factorization theorem |13] predicts that fragmentation is independent of the 
quantum numbers of both the projectile and target. Thus one expects ^ = to 
leading order. However, it is possible that the forward-moving heavy quarks will 
coalescence with the spectator valence quarks of the projectile to produce leading 
hadrons in the final state. In a gauge theory one expects the strongest attraction to 



occur when the spectator and produced quarks have equal velocities [Q. Thus the 
coalescence probability should be largest at small relative rapidity and relatively low 
transverse momentum where the invariant mass of the Qq system is minimal, and its 
amplitude for binding is maximal. 

The coalescence of charmed quarks with projectile valence quarks may also occur 
in the initial state. For example, the vr" or p wavefunctions can fluctuate into \udcc) 
or \uudcc) Fock states. These states are produced in QCD from amplitudes involving 
two or more gluons attached to the charmed quarks. The most important fluctuations 
occur at minimum invariant mass M where all the partons have approximately the 
same velocity. These fluctuations can have very long lifetimes in the target rest 
frame, (9(2Piab/M^), where Piab is the projectile momentum. Since the charm and 
valence quarks have the same rapidity in these states, the heavy quarks carry a large 
fraction of the projectile momentum. Furthermore the comoving heavy and light 
quarks can readily coalesce to produce leading charm correlations at a large combined 
longitudinal momentum. Such a mechanism can dominate the hadroproduction rate 
at large xp- This is the underlying assumption of the intrinsic charm model [0. 



The intrinsic charm fluctuations in the wavefunction are initially far off the light- 
cone energy shell shell. However, they become on shell and materialize into charmed 
hadron when a light spectator quark in the projectile Fock state interacts in the target 
T^. Since such interactions are strong, the charm production will occur primarily 



on the front face of the nucleus in the case of a nuclear target. Thus an important 
characteristic of the intrinsic charm model is its strong nuclear dependence; the cross 
section for charm production via the materialization of heavy Fock states should have 
a nuclear dependence at high energies similar to that of inelastic hadron-nucleus cross 
sections. 

In this work, we concentrate on the charmed hadrons and meson pairs channels 
studied by E791 in order to further examine the relationship between fragmentation 
and coalescence mechanisms. The calculations are made within a two-component 



model: leading-twist fusion and intrinsic charm ||T^, |T^, |18| . We find that the coales- 
cence of the intrinsic charmed quarks with the valence quarks of the projectile is the 
dominant mechanism for producing fast D and D* mesons. On the other hand, when 
the charmed quarks coalesce with sea quarks, there is no leading charmed hadron. 
We discuss the longitudinal momentum distributions and the related asymmetries for 
Ac/Ac and Dg/Dg production, as well as DD pairs. (We have only applied our model 
to DD pairs in the forward hemisphere in order to provide a clear definition of the 
asymmetry.) 

As expected, the asymmetries predicted by the intrinsic charm coalescence model 
are a strong function of xp. We find that Ac production in the proton fragmentation 
region (xi;' < in n^p collisions) is dominated by the coalescence of the intrinsic 
charm quark with the ud valence quarks of the proton. Coalescence is particularly 
important in DD pair production. The production of Dg/Dg and, at xp > 0, Ac/ Ac 
by coalescence must occur within still higher particle number Fock states. 

Leading particle correlations are also an integral part of the Monte Carlo program 
PYTHIA [|19i based on the Lund string fragmentation model. In this model it is 



assumed that the heavy quarks are produced in the initial state with relatively small 
longitudinal momentum fractions by the leading twist fusion processes. In order 
to produce a strong leading particle effect at large xp, the string has to accelerate 
the heavy quark as it fragments and forms the final-state heavy hadron. Such a 
mechanism goes well beyond the usual assumptions made in hadronization models 
and arguments based on heavy quark symmetry, since it demands that large changes 
of the heavy quark momentum take place in the final state. 

In this paper we shall compare the predictions of the intrinsic charm coalscence 
model with those of PYTHIA The comparison of the data with these models 
should distinguish the importance of higher heavy quark Fock state fluctuations in 
the initial state and the coalescence process from the strong string hadronization ef- 
fects postulated in the PYTHIA model. 



2. Leading-Twist Production 

In this section we briefly review the conventional leading-twist model for the 
production of single charmed hadrons and DD pairs in 7r~p interactions. We will also 
show the corresponding distributions of charmed hadrons predicted by the PYTHIA 
model [0 . 

Our calculations are at lowest order in a^. A constant factor ~ 2 — 3 is included 
in the fusion cross section since the next-to-leading order xp distribution is larger 
than the leading order distribution by an approximately constant factor . Neither 
leading order production nor the next-to-leading order corrections can produce flavor 
correlations 

The single charmed hadron xp distribution, xp = {2mT/ \/s) sinh?/, has the fac- 



torized form [T^ 



- -TT / Hab{Xa,Xb)— ' dzs dy2 dpp , (2) 



dxi 



where a and h are the initial partons, 1 and 2 are the charmed quarks with irtc = 1.5 
GeV, and 3 and 4 are the charmed hadrons. The convolution of the subprocess cross 
sections for qq annihilation and gluon fusion with the parton densities is included in 

HabiXa, Xb) = T.[fgi^a)ff{Xb) + (x.) /f (x,)] ^ + ff{x^)f^{x,)^ , (3) 

where A and B are the interacting hadrons. For consistency with the leading-order 
calculation, we use current leading order parton distribution functions, GRV LO, for 



both the nucleon ||22| and the pion p3|] . 

The fragmentation functions, Dh/c{z), describe the hadronization of the charmed 
quark where z = Xh/xc is the fraction of the charmed quark momentum carried 
by the charmed hadron, produced roughly collinear to the charmed quark. Assuming 
factorization, the fragmentation is independent of the initial state (leptons or hadrons) 
and thus cannot produce flavor correlations between the projectile valence quarks 
and the charmed hadrons. This uncorrelated fragmentation will be modeled by two 
extremes: a delta function, 6{z — l), and the Peterson function [^, as extracted from 
e^e~ data. The Peterson function predicts a softer xp distribution than observed in 
hadroproduction, even at moderate xf [1^, since the fragmentation decelerates the 



charmed quark, decreasing its average momentum fraction, (xp), approximately 30% 
by the production of D mesons. The delta-function model assumes that the charmed 
quark coalesces with a low-x spectator sea quark or a low momentum secondary 



quark such that the charmed quark retains its momentum This model is more 
consistent with low px charmed hadroproduction data ||2^, |2B|, ^ than Peterson 
fragmentation. 

The parameters of the Peterson function we use here are taken from e"'"e~ studies 
of D production [^. The D* distributions are very similar to the D distributions 
but the Ac distribution appears to be somewhat softer Although there is some 
uncertainty in the exact form of the Peterson function for charmed baryons and Dg 
mesons, it always produces deceleration. 
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Figure 1: Leading- twist fusion calculations of single charm production from vr p inter- 
actions at 500 GeV with delta function (a) and Peterson function (b) fragmentation. 



The results from the PYTHIA Monte Carlo ^ for AJAc (c) and (d) pro- 

duction are also shown. The fusion calculations are normalized to the charm cross 
section while the PYTHIA distributions are normalized to the number per event. 

In Fig. 1 we show the single inclusive xp distributions calculated for (a) delta 
function and (b) Peterson function fragmentation in 7T~p interactions at 500 GeV. The 
results are normalized to the total single charmed quark cross section. The parton 
distributions of the pion are harder than those of the proton at large Xbj, producing 
broader forward distributions. As expected, the delta-function fragmentation results 
in harder distributions than those predicted by Peterson fragmentation for xp > 0.2. 
However, as shown in |Q, the conventional fusion model, even with delta-function 
fragmentation, cannot account for the shape of leading D distributions. 

The charmed hadron distributions from PYTHIA, obtained from a ir^p run with 
10® events at 500 GeV using all default settings and the GRV LO parton distribu- 
tions, are shown in Fig. 1(c) for A^ and A^ and 1(d) for Dg and Dg hadrons. The 
distributions are normalized to the number of charmed hadrons per event. PYTHIA 
is based on the Lund string fragmentation model fl^ in which charmed quarks are 
at string endpoints. The strings pull the charmed quarks toward the opposite string 
endpoints, usually beam remnants. When the two string endpoints are moving in the 
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same general direction, the charmed hadron can be produced with larger longitudinal 
momentum than the charmed quark, accelerating it. In the extreme case where the 
string invariant mass is too small to allow multiple particle production, this picture 
reduces to final-state coalescence and the string endpoints determine the energy, mass, 
and fiavor content of the produced hadron Thus a D~ or can inherit all of 



the remaining projectile momentum while and -0° production is forbidden. The 
coalescence of a charmed quark with a valence diquark results in the secondary peak 
at xp ~ —0.85 in the distribution shown in Fig. 1(c). The baryon is a leading 
charmed hadron in the proton fragmentation region since it can have two valence 
quarks in common with the proton. Also in the proton fragmentation region, the Ds 
is somewhat harder than the Ds- This is evidently a secondary effect of K{uds) coa- 
lescence with a ud diquark: the s can pull the charmed quark to larger Xi? in the wake 
of the A, producing more Dg at negative xp than Dg- Such coalescence correlations 
are not predicted for the A^ or the in the proton fragmentation region. In the pion 
fragmentation region, the predicted Dg/Ds distributions are harder than the Ac/Ac 
distributions. 

The charmed pair xp distribution is 

= / Hab{Xa,Xi,)^^ ^^''^ "'"^ ^(M^;p-2m|(l + cosh(|/3-|/4))) 
dXp J il/\il/2 z^z^ 

X S{xp — X3 — X4) dz^dzidy^dyiidppdM'jj-^ , (4) 



where xp = x^ + x^ and m\ = p^ + m^). Figure 2 shows the forward xp distribution of 
DD pairs from TT~p interactions at 500 GeV with (a) delta function and (b) Peterson 
function fragmentation for both charmed quarks. The distributions are normalized to 
the total cc pair production cross section. Note that the pair distributions are harder 
than the single distributions in Fig. 1. The perturbative QCD calculation cannot 
distinguish between leading and nonleading hadrons in the pair distributions. 

In Fig. 2(c), the DD pair distributions from PYTHIA have been classified as 
doubly leading (LL), D~D^, nonleading-leading (NL), D^D^ + D^D'^, and doubly 

9 



10-^ I ' ' ' ' I ' ' ' ' 
PYTHIA (c) 

_ , ° LL °Ooooo„ □ 

-D 10--" h*NL 

°NN DD 

I I I I I I I I !□ 



10" 



I ' ' ' ' 

PYTHIA (d) 



90ooHHgnn?e 



D*D* 

. I . 







12-95 
8096A2 



0.5 1.0 

XF 



0.5 1.0 

XF 



Figure 2: Leading-twist fusion calculations of cc pair production from vr p interactions 
at 500 GeV with delta function (a) and Peterson function (b) fragmentation. The 



results from the PYTHIA Monte Carlo [IS] for LL, NL, and A^A^ DD (c) and DW 



(d) pairs are also shown. The fusion calculations are normalized to the cc cross section 
while the PYTHIA distributions are normalized to the number per event. 

nonleading (NN), D^D^ . The same classification for D* pairs is shown in Fig. 2(d). 
The distributions are normalized to the number of charmed pairs per event. We have 
not considered DD* or D*D pairs. Note that the leading particle assignments are 
only valid for xp > 0. The assignments are more meaningful for the D*'s since they 
may be assumed to be directly produced. The neutral D mesons in Fig. 2(c) arise 
in part from charged D* decays. Thus the A^A^ D pair distributions have a shoulder 
at xp ~ 0.8 from D*^ decays to D'^ which is absent in the A^A^ D*D* distributions. 
Note also that the NL pairs are most numerous since both charged and neutral D 
pairs contribute to the distribution. Almost three times as many neutral D's are 
produced than charged D's. Only some of this difference can arise from D* decays 
since charged and neutral D*^s are produced in nearly equal abundance. 

Other final-state coalescence models have been proposed, including a valence spec- 
tator recombination model 0] and the valon model Two important unknowns 



in these models are the correlation between the charmed quark and the valence spec- 
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tator in the recombination function and the ra-particle parton distributions of the 
spectator and participant valence quarks. In this work we will not compare to either 
of these models but simply note that they can also produce charmed hadrons and 
hadron pairs by final-state coalescence. 



3. Intrinsic Heavy Quark Production 

The wavefunction of a hadron in QCD can be represented as a superposition of 
Fock state fluctuations, e.g. \nv), \nvg), \nvQQ), ...components where ny = ud 
for a TT~ and uud for a proton. When the projectile scatters in the target, the co- 
herence of the Fock components is broken and the fluctuations can hadronize either 
by uncorrelated fragmentation or coalescence with spectator quarks in the wavefunc- 



tion [15, [T^. The intrinsic heavy quark Fock components are generated by virtual 
interactions such as gg QQ where the gluons couple to two or more projectile 
valence quarks. The probability to produce QQ fluctuations scales as a'j.{mQQ)/mQ 
relative to leading- twist production and is thus higher twist. Intrinsic QQ Fock 
components are dominated by configurations with equal rapidity constituents so that, 
unlike sea quarks generated from a single parton, the intrinsic heavy quarks carry a 
large fraction of the parent momentum [l^ . 



The frame-independent probability distribution of an n-particle cc Fock state is 



dxi--- dxn ' '''' {m\ - Y^'l=i(rh'f/xi)y ' 

where Nn, assumed to be slowly varying, normalizes the \nvcc) probability. Pic. The 
delta function conserves longitudinal momentum. The dominant Fock configurations 
are closest to the light-cone energy shell shell and therefore have minimal invariant 
mass, = J^i^l/^iy where m? = k^^ + mf is the effective transverse mass of the i^^ 
particle and Xi is the light-cone momentum fraction. Assuming (fc^j) is proportional 
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to the square of the constituent quark mass, we adopt the effective values rriq = 0.45 
GeV, m, = 0.71 GeV, and m, = 1.8 GeV ITZl, HI. 



The intrinsic charm production cross section can be related to Pic and the inelastic 
hN cross section by 

aUhN) = P-A^^ . (6) 

The factor of /i^/4m^ arises because a soft interaction is needed to break the coherence 
of the Fock state. The NA3 collaboration |]3^ separated the nuclear dependence of 



J/ip production in tt~ A interactions into a "hard" contribution with a nearly linear 
A dependence at low xp and a high xp "diffractive" contribution scaling as A^ ''^, 
characteristic of soft interactions. One can fix the soft interaction scale parameter, 
~ 0.2 GeV^, by the assumption that the diffractive fraction of the total production 
cross section is the same for charmonium and charmed hadrons. Therefore, we 
obtain cric(7rA^) ^ 0.5 /ib at 200 GeV and (T;^{pN) ^ 0.7 /ib [|12| with Pjc = 0.3% from 



an analysis of the EMC charm structure function data ^5]. A recent reanalysis of 
this data with next-to-leading order calculations of leading twist and intrinsic charm 
electroproduction confirms the presence of an 1% intrinsic charm component in 



the proton for large xbj |^6|. Note that a larger Pic would not necessarily lead to a 
larger a-^c- Since we have fixed /i^ from the NA3 data, increasing Pic would decrease 
/i^ accordingly. 

We now calculate the probability distributions, dPic/dxp for charmed hadrons 
and DD pairs resulting from both uncorrelated fragmentation and coalescence of the 
quarks in the intrinsic charmed Fock states. These light-cone distributions are frame 
independent. In a hadronic interaction, these states are dissociated and materialize 
with the corresponding differential cross section 

dxF ~ ""^^AmldxF ■ ^ ^ 

In the case of vr^p collisions, the fluctuations of the \udcc) state produces charmed 
hadrons at Xi;' > in the center of mass while the fluctuations of the \uudcc) state 
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produces charmed hadrons ai xp < 0- 

3.1 Single charmed hadrons 

There are two ways of producing charmed hadrons from intrinsic cc states. The 
first is by uncorrelated fragmentation, discussed in Section 2. Additionally, if the pro- 
jectile has the corresponding valence quarks, the charmed quark can also hadronize by 
coalescence with the valence spectators. The coalescence mechanism thus introduces 
flavor correlations between the projectile and the final-state hadrons, producing e.g. 
D~'s with a large fraction of the tt" momentum. In the pion fragmentation region, 
xp > 0, D~ and have contributions from both coalescence and fragmentation 
while and can only be produced from the minimal cc Fock state by fragmen- 
tation. In the proton fragmentation region, xp < 0, the D~ and D° may be produced 
by both coalescence and fragmentation. 

If we assume that the c quark fragments into a D meson, the D distribution is 



where n — 4, 5 for pion and proton projectiles in the \nvcc) configuration. This 
mechanism produces D mesons carrying 25-30% of the projectile momentum with 
the delta function and 17-20% with the Peterson function. The D distributions are 
shown in Fig. 3(a) and 3(b) for proton and pion projectiles normalized to the total 
probability of the \nvcc) Fock state configuration assuming that Pic is the same for 
protons and pions. Less momentum is given to the charmed quarks in the proton 
than in the pion because the total momentum is distributed among more partons. 
These distributions are assumed for all intrinsic charm production by uncorrelated 
fragmentation. 

The coalescence distributions, on the other hand, are specific for the individual 
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Figure 3: Charmed particle distributions from the intrinsic charm modeL The 
charmed quark fragmentation distributions are shown for Peterson function (sohd) 
and delta function (dashed) fragmentation for proton (a) and pion (b) projectiles. 
Charmed hadrons produced by coalescence from protons (solid) and pions (dashed) 
are also shown for (c) Ac baryons and (d) Ds mesons. The charmed quark distribu- 
tions in (a) and (b) and the distribution from a proton (c) are normalized to Pic. 
The Ac distribution from a pion (c) and the Dg distributions (d) are normalized to 
Picu and Pics respectively. 
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charmed hadrons. The coalescence contribution to leading D production is 

-P^ = f \\dxi ^^^^ — 5{xd-Xc-xi). (9) 
dxu J j^J[ axi . . . axn 

With the additional momentum of the light valence quark, the D takes 40-50% of 

the momentum. In the proton fragmentation region, the c quark can coalesce with 

valence u and d quarks to produce leading A^'s, 

■IT^ = lid^i^JZ '^SixA.-^c-Xi-X2) , (10) 

carrying 60% of the proton momentum. The distribution, shown in Fig. 3(c), is also 
normalized to Pic. 

Coalescence may also occur within higher fluctuations of the intrinsic charm Fock 
state. For example, at > 0, Ac and Dg can be produced by coalescence from 
\nvccdd) and \nvccss) configurations. We previously studied ipip production from 
\nvcccc) states [^. Assuming that all the measured ipip pairs ||3^, |3^ arise from 



these configurations, we can relate the ipip cross section. 



4''{hN) = fi/,^aUhN), (11) 



IC 



to the double intrinsic charm production probability. Pice, where f^/h is the fraction 
of intrinsic cc pairs that become J/ip^s. The upper bound on the model, a^^ = 



iV) ^ 20 pb m, requires Pice ~ 4.4% Pic |0[. This value of Pice can be 



used to estimate the probability of light quark pairs in an intrinsic charm state. We 
expect that the probability of additional light quark pairs in the Fock states to be 
larger than Pice, 

leading to Pjeu = Pied ~ 70.4% Pje and Pjes ~ 28.5% Pje. 

Then with a pion projectile at xp > 0, the Ac coalescence distribution from a 
six-particle Fock state is 

= / n ^^^T~~^ — ^^^^'^ -Xc-xi- X2) , (13) 
15 



also shown in Fig. 3(c) and normalized to Picu- Since half of the quarks are needed to 
produce the it carries 50% of the pion momentum. The Dg mesons arising from 

coalescence in the \nvccss) state, 

carry ^ 30-40% of the hadron momentum, as shown in Fig. 3(d) and normalized to 
Pics- The Dg and Dg inherit less total momentum than the leading D since the Fock 
state momentum is distributed over more partons. Thus as more partons are included 
in the Fock state, the coalescence distributions soften and approach the fragmentation 
distributions, Eq. (8), eventually producing charmed hadrons with less momentum 
than uncorrelated fragmentation from the minimal cc state if a sufficient number of 
qq pairs are included. There is then no longer any advantage to introducing more 
light quark pairs into the configuration-the relative probability will decrease while 
the potential gain in momentum is not significant. We thus do not consider Ac pro- 
duction by coalescence at xp < since a minimal nine-parton Fock state is required. 



3.2 DD pair production 

In any |nycc) state, DD pairs may be produced by double fragmentation, a com- 
bination of fragmentation and coalescence, or, from a pion projectile only, double 
coalescence. We discuss only pair production at > so that \udcc) is the minimal 
Fock state. In the proton fragmentation region, with no valence antiquark, the lead- 
ing mesons would be D'(cd) and D^{cu). Therefore no doubly leading D~D^ pairs 
can be produced from the five parton state: two intrinsic cc pairs are needed, i.e. 
\uudcccc) states, automatically softening the effect. However, doubly leading meson- 
baryon pairs such as AcD^ and E^+D" may be produced by coalescence in the \uudcc) 
state. These combinations might be interesting to study in pp interactions. 
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Figure 4: Intrinsic DD distributions from a vr^ projectile. Double fragmentation 
distributions are shown for Peterson function (solid) and delta function (dashed) 
fragmentation in (a). Charmed pairs produced by coalescence of one quark and 
fragmentation of the other are shown in (b) for Peterson (solid) and delta (dashed) 
function fragmentation. The double coalescence production from a pion is a delta 
function at = 1. The distributions are normalized to Pic. 
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The DD pairs resulting from double fragmentation, 
dPi^^ /• , , ^ , dP 



dxjjj-) 



/dr\c 
dZcdZc n ^^ Vxi '^dXn ^^'"^^"^^^''^^^^^ ^^^^ 
X 5{xd - ZcXc)5{x^ - ZcXc)S{Xjj-^ - xd - x^) , 



carry the lowest fraction of the projectile momentum. The distributions are shown in 
Fig. 4(a). When both charmed quarks fragment by the Peterson function, ~ 40% of 

the TT~ momentum is given to the pair while the average is 62% with delta function 
fragmentation. If, e.g. a D is produced by coalescence while the D is produced by 
fragmentation, 

= dzlldxi " D^/-{z) (16) 

^•^DD 1=1 dXi . . . UXji 

X 5{xd - Xc- Xi)5{x^ - ZXc)S{Xjy^ - Xd - x^) . 

These distributions, with rather large momentum fractions, 71% for Peterson frag- 
mentation and 87% for the delta function, are shown in Fig. 4(b). When the projectile 
has a valence antiquark, as in the pion, DD pair production by double coalescence is 
possible, 

dxjjf^ J dxi . . . dxji 

All of the momentum of the four-particle Fock state is transferred to the DD pair, 
i.e. Xjj^ = 1. We also consider double coalescence from a pion in a six particle Fock 
state. In this case, 74% of the pion momentum is given to the pair, similar to the 
result for Peterson fragmentation with coalescence, Eq. (16), as could be expected 
from our discussion of Ds production in this model. 



/ n ^^'dx '^dx~^^^^ - Xc- xi)5{x^ -Xc- X2)5{xjj-^ -Xd- x^) .(17) 



4. Predictions of the Two- Component Model 

We now turn to specific predictions of the xp distributions and asymmetries in 
our model. The Xf distribution is the sum of the leading-twist fusion and intrinsic 



charm components, 



da da\t da\c 



dxp dxp dxp 

where daiddxp is related to dP^c/dxp in Eq. (7). Note that when we discuss uncorre- 
lated fragmentation, the same function, either the deha or Peterson function, is used 
for both leading twist fusion and intrinsic charm. The intrinsic charm model pro- 
duces charmed hadrons by a mixture of uncorrelated fragmentation and coalescence 
| ]T^ 0. Coalescence in the intrinsic charm model is taken to enhance the leading 
charm probability over nonleading charm. Since we have not made any assumptions 
about how the charmed quarks are distributed into the final-state charmed hadron 
channels, an enhancement by coalescence is not excluded as long as the total probabil- 
ity of all charmed hadron production by intrinsic charm does not exceed Pic- Because 
little experimental guidance is available to help us separate the charm production 
channels, we have not directly addressed the issue here. Thus the above distributions 
are normalized to the total cc cross section for the pair distributions and to the single 
charm cross section for the single charmed hadrons. This is naturally an overestimate 
of the cross sections in the charm channels, but more complete measurements are 
needed before the relative strengths of the D, A^, Dg, D*, etc. contributions to the 
charm cross section can be understood. 

In the case of nuclear targets, the model assumes a linear A dependence for leading- 
twist fusion and an ^4° '''' dependence for the intrinsic charm component This A 
dependence is included in the calculations of the production asymmetries while the 
xp distributions are given for 7c~p interactions. The intrinsic charm contribution to 
the longitudinal momentum distributions is softened if the A dependence is included. 



4.1 Single charmed hadrons 

We now consider the single charmed hadron distributions produced in 7r~p in- 
teractions at 500 GeV over the entire xp range. Since the production mechanisms 
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are somewhat different for positive and negative xp, particularly for the Ac, we will 
discuss the pion and proton fragmentation regions separately. 

We begin with Ac production in the proton fragmentation region, negative xp. As 
we have already indicated, the Ac can be produced by coalescence from the \uudcc) 
configuration with an average of 60% of the center-of-mass proton momentum. The Ac 
can only be produced by fragmentation from a five-particle Fock state and, if a nine- 
particle Fock state is considered, the coalescence distribution will not be significantly 
harder than the fragmentation distribution shown in Fig. 3(a) since four additional 
light quarks are included in the minimal proton Fock state. 

Therefore coalescence is only important for the Ac, leading naturally to an asym- 
metry between Ac and Ac. We will assume that the same number of Ac and Ac are 
produced by fragmentation and that any excess of Ac production is solely due to 
coalescence. Then, at < 0, 

dPt 

dPt 

dxK^ 

The fragmentation distribution is taken from Eq. (8), the coalescence distribution 
from Eq. (10). The parameter r is related to the integrated ratio of Ac to Ac pro- 
duction. We have assumed three values of r: 1, 10, and, as an extreme case, 100. 
The results for the two uncorrelated fragmentation functions are shown in Fig. 5(a) 
and 5(b). Intrinsic charm fragmentation produces a slight broadening of the Ac dis- 
tribution for delta function fragmentation over leading-twist fusion (increasing to a 
shoulder for the Peterson function). The Ac distribution, strongly dependent on r, is 
considerably broadened. 

The value r = 1 is compatible with early low statistics measurements of charmed 
baryon production where equal numbers of Ac and Ac were found in the range 
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Figure 5: The A^/Ac xp distributions predicted by the two-component model in (a) 
and (b) for delta and Peterson function fragmentation respectively. The associated 
asymmetry is shown in (c). The xp distributions are normalized to our calculated 
cross section, Eq. (18). In (a) and (b), the solid curve is the Ac distribution (identical 
to Ac for xp > 0) while the dashed, dot-dashed, and dotted curves show Ac distribu- 
tions with r = 1, 10, and 100. The predicted model asymmetry for r = 1 (solid), 10 
(dashed), and 100 (dot-dashed) is compared with that from PYTHIA (open circles) 
in (c). At xp > 0, our model predicts no asymmetry. 
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1x^1 < 0.3. The data is often parameterized as (1 — where 

1 — Pi? mill „ /oi N 

= n ^ r - 2 . (21) 



For xp < 0, we predict = 4.6 for the delta function and 7.3 for the Peterson 
function. The difference is due to the steeper slope of the fusion model with the 
Peterson function. We find rather large values of since the average Xi;' is dominated 
by the leading- twist fusion component at low xp- If we restrict the integration to 
xp < —0.5, then decreases to 1.4 independent of the fragmentation mechanism^. 
The data on Ac production measured in pp collisions at the ISR with ^/s = 63 GeV 
are consistent with this prediction. For xp > 0.5, n\^, = 2.1 ± 0.3 was found while 
for xp > 0.35, riAc = 2.4 ± 1.3 0. Hard charmed baryon distributions have also 
been observed at large xp in nN interactions at the Serpukhov spectrometer with 
an average neutron energy of 70 GeV, = 1.5 ± 0.5 for xp > 0.5 [H. Charmed 



hyperons Hc('usc) produced by a 640 GeV neutron beam do not exhibit a strong 
leading behavior, ris^ = 4.7 ± 2.3. This is similar to the delta function prediction for 
when xp < 0. On the other hand, charmed hyperons produced with a T,~{dds) 
beam |TD|, |rT| are leading with n^^ = 1.7 ± 0.7 for xp > 0.6 [|T^]. Thus in the 
proton fragmentation region r = 1 is compatible with the shape of the previously 
measured Ac xp distributions. When we compare the Ac cross section in the proton 
fragmentation region with that of leading-twist fusion, the coalescence mechanism 
increases the cross section by a factor of 1.4-1.7 over the fusion cross section and by 
30% over the Ac cross section. 

The extreme value, r = 100, was chosen to fit the forward Ac production cross 
section measured at the ISR, Bapp^\^x = 2.84±0.5 0, assuming that the charmed 
quark and Ac cross sections are equal, already an obvious overestimate. This choice 
produces a secondary peak in the 7T~p distributions at xp ~ —0.6, the average Ac 



^The parameterization (1 — xp)^ is only good if the distribution is monotonic. However, our 
two-component Ac distribution does not fit this parameterization over all xp. At low xp, the 
leading-twist component dominates. If only the high xp part is included, the value of ([xfI) is a 
more accurate reflection of the shape of the intrinsic charm component. 



22 



momentum from coalescence. Such a large value of r implies that the intrinsic charm 
cross section is considerably larger than the leading-twist cross section. 

The shape of the distribution with r = 100 is similar to that due to diquark 



coalescence in PYTHIA |T^, shown in Fig. 1(c), except that the PYTHIA distribu- 
tion peaks at xp ~ —0.9 due to the acceleration induced by the string mechanism. 
While a measurement of the Ac cross section over the full phase space in the proton 
fragmentation region is lacking, especially for pp interactions at Xi? > 0, no previous 
measurement shows an increase in the A^ x^? distributions as implied by these results. 
However, the reported Ac production cross sections are relatively large |^, ^, 
between 40 /xb and 1 mb for 10 < < 63 GeV. In particular, the low energy cross 
sections are much larger than those reported for the cc total cross section at the same 
energy. This is not yet understood. 

A few remarks are in order here. Some of these experiments [0, H extract the total 
cross section by extrapolating flat forward xp distributions back to xp = ^ and also 
assume associated production, requiring a model of D production. On the other hand, 
the reported cc total cross sections are usually extracted from D measurements at 
low to moderate xp and would therefore hide any important coalescence contribution 
to charmed baryon production at large xp. High statistics measurements of charmed 
mesons and baryons over the full forward phase space {xp > 0) in pp interactions 
would help resolve both the importance of coalescence and the magnitude of the total 
cc production cross section. 

We also chose r = 10 as an intermediate value. In this secondary peak is 

also predicted but the cross section at x^? < is only a factor of two to three larger 
than the fusion cross section rather than the factor of 21 needed to fit the ISR data 
at Xp > 0.5. The magnitude of the second maximum is also less than the fusion cross 
section in the central region. 

An important test of the production mechanism for charm hadroproduction is the 
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Ac and Ac asymmetry, defined as 

_ _ da{K) - da(K,) 
^^^/^^ rfa(Ac) + rfa(Ac) ■ ^ ^ 

If A^^ij^^ is assumed to arise only from initial state coalescence, we can estimate the 
parameter r from the E791 500 GeV t[~ A data. Our calculated asymmetries for the 
three r value^ are compared with the results from PYTHIA in Fig. 5(c). It is, as 
expected, closest to our model with r = 100 although the asymmetry predicted by 
PYTHIA does not increase as abruptly as in our model. 

Preliminary data from E791 indicate a significant asymmetry for Xi? as small 



as —0.1, albeit with large uncertainty. The intrinsic charm model in its simplest 
form can only produce such asymmetries if r > 100, against intuition. Alternatively, 
a softer Ac distribution from coalescence would make a larger asymmetry at lower 
\xp\, thus allowing a smaller r. Such a softening could be due to either an important 
contribution to Ac production from a \nvgcc) configuration or a different assumption 
about the \nycc) wavefunction |Q. In a more realistic model, both initial and final 
state coalescence will play some role in Ac production. Final-state coalescence, added 
to the leading-twist fusion prediction, would also require a smaller r. 

At Xi? > there is no asymmetry in 7r~p interactions since both the baryon and 
antibaryon can be produced by fragmentation from a \udcc) state and by coalescence 
from a \udccqq) state (g = u^d). Then 

dP-^" dP^" dP-^ P dP^ 

ic _ ic _ ic _|_ ^ icq ic ^23^ 

dxK^ dxji^ dxA^ Pic dxAc 
The coalescence contribution, obtained from Eq. (13), produces a small shoulder in 
the distributions at xp > 0. We extract = 4.1 for xpmin = 0, in good agreement 



with the NA32 measurement, 12^^ = 3.5 ±0.5 ||2^. The same mechanism can account 
for both Ac and Ac production in the tt~ fragmentation region since no asymmetry 
is observed [43|, which is also in accord with the NA32 result, cr(Ac)/cr(Ac) ~ 1 |27|. 



"''We have only shown the delta function results. Those with the Peterson function are quite 
similar. The slope increases slightly but the point where /a > docs not shift. 
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To look for subtle coalescence effects as well as to understand the large difference 
between and Ac production at xp ~ —0.1 it is important to compare the shape of 
the momentum distributions in addition to the asymmetries. 

Assuming that equal numbers of Dg and Dg are produced, the xp distributions 

are 

dP^" dP^' dP-^ P dP^ 

ic _ ic _ IC _|_ ICS ic /2^N 

dxD, dx^^ dxD, Pic dxD, 
The coalescence contribution is given by Eq. (14) and the distributions are shown in 
Fig. 6. The shoulder at a;^ > predicted in charmed baryon production is absent 
for Ds production. Coalescence does not produce a significant enhancement of Dg's 
since Pics < -Picu- The average momentum gain over uncorrelated fragmentation of 
the \nvcc) state is small. The forward distributions in Fig. 6 are only slightly harder 
than those from leading-twist production. This is also true for xp < where the 
Dg, Dg distributions are not significantly different from the Ac distributions even 
though coalescence is included in the production of the Dg and not in Ac. When we 
compare our distribution with the parameterization (1 — Xi7')"^% we extract ud^ = 4.7 



at xp > 0, in agreement with the NA32 measurement, = 3.9 ± 0.9 [27 



We find -^^s/Ds (^^) ~ ^ -^^ model since the production mechanisms 

are everywhere identical for the particle and antiparticle. In contrast to the intrinsic 
charm model, the Dg excess predicted by PYTHIA in the proton region, shown in 
Fig. 1(d), leads to a backward asymmetry similar to A^^^j^^, shown in Fig. 6(c). 



4.2 DD production 

We simplify our discussion of DD pair production by several respects. We have 
assumed that equal numbers of D* mesons and, separately, primary D mesons are 
produced by fragmentation and that any production excess is the result of coalescence. 
Thus the DD and D*D* xp distributions are equivalent within the model and the 
analysis applies for both types of pairs unless stated otherwise. Therefore we shall 

25 




Figure 6: The DgjDs xp distributions predicted by the two- component model in (a) 
and (b) for delta and Peterson function fragmentation respectively. The associated 
asymmetry from our model and PYTHIA (open circles) is shown in (c). Our model 
predicts no asymmetry. 
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also implicitly assume that all secondary D^s produced by D* decays can be separated 
from the primary D's. We do not consider DD* or D*D pairs. 

We use the same pair classification as in our discussion of the pair distributions 
from PYTHIA, shown in Fig. 2(c) and 2(d). Then the probability distributions for 
intrinsic charm pair production are: 



NN A tdFF 

(25) 



dxj-)j-) dxj^j-, 



dP-^^ fdP-^^ dP^^'\ 

p£_ = 2 p^ + 1.2p2_ (26) 

d^DD \ d^DD dx^^ I 

jpLL jpFF fdP'"^ P dP'^^\ 

p3 = + 5^ + P,^^5(x^5-l) + ^i^ . (27) 

'^•^ D D "'•^ DD \ '^■^ D D -'ic DD J 

In the above, we have assumed that there is a 20% production enhancement due to 
coalescence, e.g. 



dPP~ dP-^ dP'^ 

= ^ + 1.2^ . (28) 
dxn- dxn- dxo- 



This assumption is different from our calculation of Ad-/d+ in IHlj where we assumed 
that the D~ and D^ production probabilities were equal even though the D~ is pro- 
duced by coalescence and the D^ is not. This resulted in a small negative asymmetry 
at low xp. When Ad-/d+ is recalculated with the D~ distribution in Eq. (28), the 
model asymmetry is never negative and better agreement with the data [0, ^, Q is 
found. 

A factor of two has been included in the NL distribution because there are two 
sources of the NL pairs {D^D^ and D^D^) relative to the LL and A^A^ pairs. This 
factor is also included in the leading-twist production cross section. The last term 
in the LL distribution, from the six-particle Fock state, is the only significant source 
of pairs due to double coalescence, causing the LL distributions to be somewhat 
harder than the NL distributions. The pair distributions are given in Fig. 7(a) and 
7(b). Note that even the A^A^ distribution with the Peterson function in Fig. 7(b) is 
broadened considerably over the leading twist distribution shown in Fig. 2(b). The 
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predicted slopes of the pair distributions from the various sources begin to differ for 
xf > 0.25. 




Figure 7: The DD pair distributions from 7c~p interactions at 500 GeV for A^A^ 
(sohd), NL (dashed), and LL (dot-dashed) pairs are shown in (a) and (b) for delta 
function and Peterson function fragmentation. The intrinsic charm model distribu- 
tions are given in Eqs. (25)- (27). The asymmetries for delta function and Peterson 
function fragmentation from the two- component model, calculated using Eqs. (29)- 
(31), are shown in (c) and (d). They are compared with the D*D asymmetries from 
PYTHIA for All-nn (solid curve and squares), All-nl (dashed curve and circles), 
and Anl-nn (dot-dashed curve and crosses). 

We define the following three asymmetries: 

dajLL) - dajNN) 
•^^^-^^ - daiLL) + daiNN) ^^^^ 

dajLL) - dajNL) 
-^^^-^^ - daiLL) + daiNL) 
_ dajNL) - dajNN) 
^^"-^^ - daiNL) + da{NN) ' 

shown in Fig. 7(c) and 7(d) for delta and Peterson function fragmentation. The 
corresponding D*D* asymmetries from PYTHIA are also shown since the leading as- 
signments are unaffected by decays. It is interesting to study all three asymmetries. 
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Since twice as many NL pairs are produced by definition, All-nl is negative for all 
xp. The change in All-nl xp •-^ 1 in Fig. 7(d) is due to double coalescence from 
the \udcc) state. Larger values of Anl-nn and All-nn ai'c predicted by PYTHIA 
than by our model due to our different assumptions about particle production. How- 
ever, the general trends are quite similar for the final-state coalescence mechanism of 
PYTHIA and the initial-state coalescence of our model. This is perhaps not surpris- 
ing since the single D asymmetries from the two models would be quite similar if the 
same assumptions were made about the initial production ratios at a;^ = 0. 



5. Conclusions 

We have studied single charmed hadron and charmed meson pair longitudinal mo- 
mentum distributions and the related asymmetries within the intrinsic charm model 
and PYTHIA. In conventional leading- twist perturbative QCD, there is no asymme- 
try between charmed and anticharmed hadrons. In the intrinsic charm model as well 
as in PYTHIA, the asymmetry is clearly due to coalescence. However, we find that 
the asymmetries alone cannot tell the full story, especially aX xp > 0. The individual 
Xp distributions are needed over all xp to unravel the production properties of the 
charmed hadrons since only the shapes of these distributions can reveal deviations 
from the fusion predictions. 

The DD pair asymmetries are also quite interesting, particularly at high pair 
momentum. We have not considered their production at negative Xp m this paper 
due to the ambiguity in leading particle assignments. However, since only D~ or 
can be produced by coalescence from the five-particle Fock state, a study of Ac-D° or 
TicD~ pairs could prove more enlightening. 

Given our predictions aX xp < {) m 7r~p interactions, it would be quite inter- 
esting if high statistics measurements of Ac production can be made in pp collisions 
at \xp\ > to clarify charmed baryon production. Since we would predict similar 
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behavior for bottom hadron production, such studies would also be of interest. Al- 
though a rise in da jdxp with xp seems counterintuitive, such distributions have been 
observed in diffractive and production in S~Be interactions |^. Studies of 



charmed baryon production by a hyperon beam within the context of this model are 
underway. 
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